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ABSTRACT 


Experimental  studies  demonstrate  a  significant  effect  of  specimen  size,  alW  ratio 
and  prior  ductile  tearing  on  cleavage  fracture  toughness  values  (Jc)  measured  in  the 
ductile-to-brittle  transition  region  of  ferritic  materials.  In  the  lower-transition  re¬ 
gion,  cleavage  fracture  often  occurs  under  conditions  of  large-scale  3delding  but 
without  prior  ductile  crack  extension.  The  increased  toughness  develops  when  plas¬ 
tic  zones  formed  at  the  crack  tip  interact  with  nearby  specimen  surfaces  which  re¬ 
laxes  crack-tip  constraint  (stress  triaxiality).  In  the  mid-to-upper  transition  region, 
small  amounts  of  ductile  crack  extension  (often  <  1-2  mm)  routinely  precede  ter¬ 
mination  of  the  J-Aa  curve  by  brittle  fracture.  Large-scale  yielding,  coupled  with 
small  amounts  of  ductile  tearing,  magnifies  the  impact  of  small  variations  in  micro¬ 
scale  material  properties  on  the  macroscopic  fracture  toughness  which  contributes  to 
the  large  amount  scatter  observed  in  measured  Jc-values. 

Previous  work  by  the  authors  described  a  micromechanics  fracture  model  to  cor¬ 
rect  measured  Jc -values  for  the  mechanistic  effects  of  large-scale  yielding.  This  new 
work  extends  the  model  to  also  include  the  influence  of  ductile  crack  extension  prior 
to  cleavage.  Ductile  crack  extensions  of  10-15  x  the  crack-tip  opening  displacement 
at  initiation  are  considered  in  plane-strain,  finite  element  computations.  The  finite 
element  results  demonstrate  a  significant  elevation  in  crack-tip  constraint  due  to 
macroscopic  “sharpening”  of  the  extending  tip  relative  to  the  blunt  tip  at  initiation  of 
growth.  However,  this  effect  is  offset  partially  by  the  additional  plastic  deformation 
associated  with  the  increased  applied  J  required  to  grow  the  crack.  The  initial  alW 
ratio,  tearing  modulus,  strain  hardening  exponent  and  specimen  size  interact  in  a 
complex  manner  to  define  the  evolving  near-tip  conditions  for  cleavage  fractiu-e.  The 
paper  explores  development  of  the  new  model,  provides  necesseuy  graphs  and  proce¬ 
dures  for  its  application  and  demonstrates  the  effects  of  the  model  on  fracture  data 
sets  for  two  pressure  vessel  steels  (A533B  and  A515). 
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1.  Introduction 


Experimental  studies  consistently  reveal  large  scatter  in  the  measured  values  of 
cleavage  fracture  toughness  for  ferritic  steels  tested  in  the  ductile-to-brittle  temper¬ 
ature  transition  region  [see  for  example:  11, 12,  17,  27,  29,  33-35\.  Moreover,  data 
sets  of  seemingly  identical  specimens  frequently  show  a  mix  of  cleavage  fractures 
with  and  without  small  amounts  of  prior  stable  tearing;  fracture  toughness  values 
for  specimens  with  prior  growth  often  are  significantly  larger  than  toughness  values 
for  specimens  without  growth.  The  large  amount  of  scatter,  coupled  with  specimens 
that  exhibit  tearing  prior  to  cleavage  and  yet  others  that  show  no  tearing,  greatly 
complicates  the  interpretation  of  toughness  data  to  define  meaningful  values  for  ap¬ 
plication  in  fracture  assessments  of  structures. 

Scatter  in  fracture  toughness  values  (Jc)  observed  in  the  transition  region  arises 
from  the  complex  interplay  between  the  mechanistic  effects  of  crack-tip  constraint 
(stress  triaxiality)  and  the  stochastic  variability  of  cleavage  resistance  at  the  micro- 
structural  level.  Once  the  conditions  of  small-scale  yielding  no  longer  apply,  interac¬ 
tion  of  crack-tip  plastic  zones  with  overall  plastic  deformation  of  the  test  specimen 
creates  a  highly  nonlinear  relationship  between  J  and  near-tip  stress  triaxiality.  For 
specimens  experiencing  significant  levels  of  plastic  deformation,  large  increases  of  J 
must  be  imposed  to  produce  increases  in  the  near-tip  stress  fields.  This  phenomenon 
magnifies  the  impact  of  small  variations  in  the  microstructural  cleavage  resistance, 
and  contributes  a  large  measure  of  the  scatter  observed  in  fracture  data.  Dodds  and 
Anderson  [3,  7]  proposed  a  stressed-volume,  micromechanics  model  for  cleavage 
fracture  to  correct  measured  Jc-values  for  the  mechanistic  effects  of  large-scale 
Snelding.  The  model  defines  a  constraint  “corrected”  Jo-value,  denoted  JssY  in  earli¬ 
er  work,  for  each  measured  Jg-value  at  fracture.  The  corrected  Jo-value  generates 
the  same  stressed-volume  of  material  at  the  crack  tip  in  small-scale  yielding  (SSY 
with  T=0)  that  generates  at  the  crack  tip  in  the  test  specimen  at  fracture.  Under 
increased  loading  in  test  specimens,  the  ratio  of  applied  J  to  Jq  (J/ Jo)  exceeds  unity 
as  large-scale  yielding  causes  constraint  to  fall  relative  to  SSY  reference  condition. 
For  example,  the  J/  Jo  ratio  grows  much  larger  than  4-5  in  shallow  notched  bend 
specimens  (a/W<0.2). 

When  small  amounts  of  ductile  tearing  precede  failure  by  cleavage,  this  picture 
becomes  more  complex.  Stress  fields  ahead  of  a  growing  crack  differ  from  those  of  a 
stationary  crack  and  growth  exposes  additional  volumes  of  material  to  high  stresses 
as  the  tip  moves  forward.  Wallin  [34]  proposed  simple  extensions  of  his  WST  statisti¬ 
cal  model  [33]  to  include  the  first  order  effects  of  stress  field  changes  and  additional 
volume.  Recent  work  by  Varias  and  Shih  [31]  on  steady-state  crack  growth  in  SSY 
and  by  O’Dowd,  et  al.  [21]  for  transient  growth  in  SSY  demonstrates  clearly  the  de¬ 
velopment  of  high  constraint  conditions  ahead  of  the  growing  crack,  even  when 
growth  begins  under  low  constraint  conditions  (T/ao<0).  These  finite  element  anal¬ 
yses  reveal  a  significant  elevation  of  stress  triaxiality  caused  by  macroscopic  sharp¬ 
ening  of  the  extending  tip  and  by  additional  elastic  strains  as  the  crack  extends  for¬ 
ward  into  previously  unyielded  material.  In  fracture  test  specimens,  however,  these 
effects  may  be  partially  or  completely  offset  by  the  large  amount  of  plastic  deforma¬ 
tion  present  at  J/c  and  the  additional  plastic  deformation  associated  with  the  in¬ 
creased  applied  J  required  to  grow  the  crack.  The  initial  alW  ratio,  the  tearing 
modulus,  strain  hardening  exponent  and  specimen  size  interact  in  a  complex  man¬ 
ner  to  define  the  evolving  near-tip  conditions  for  cleavage  fracture. 
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This  work  extends  the  stressed-volume  constraint  model  for  cleavage  fracture  to 
include  the  mechanistic  effects  of  small  amounts  of  ductile  tearing  on  the  near— tip 
stress  fields.  Plane-strain,  finite  element  analyses  are  described  for  SE(B)  speci¬ 
mens  having  both  shallow  and  deep  cracks  (a/W  =  0.1,  0.5)  for  materials  with  hard¬ 
ening  exponents  n  =  5, 10.  The  models  are  loaded  to  prescribed  levels  of  J  =  Jjc  after 
which  crack  growth  occurs  consistent  with  the  specified  resistance  curve,  Jr  —  Act.  A 
parameter  study  for  a  range  of  Jjc  and  tearing  modulus  values  provides  the  stress 
fields  ahead  of  the  growing  cracks  needed  to  construct  J/  Jq  ratios  of  the  constraint 
correction  model.  In  this  process,  we  assume  that  ductile  crack  growth  prior  to  cleav¬ 
age  fracture  introduces  additional  complexity  only  in  the  resolution  of  near— tip 
stresses,  i.e.,  that  the  microstructural  resistance  to  cleavage  initiation  reniains 
stress-controlled  and  unaffected  by  the  small  amounts  of  prior  stable  tearing  (usual¬ 
ly  1-2  mm  or  less). 

Given  a  measured  toughness  value  at  cleavage  (Jc)  and  estimated  values  for  the 
tearing  modulus  and  Jr,  the  new  model  corrects  for  the  combined  effects  large-scale 
yielding  and  ductile  crack  extension.  Necessary  graphs  and  procedures  are  provided 
in  the  last  section  to  apply  the  model  to  experimental  data  sets.  Statistical  analyses, 
e.g.,  Weibull  methods,  developed  for  interpretation  of  cleavage  fracture  data  should 
be  applied  to  data  sets  of  corrected  e7o“Values,  with  no  distinction  then  needed  be¬ 
tween  corrected  values  for  specimens  with/without  prior  tearing.  The  paper  con¬ 
cludes  with  example  applications  of  the  proposed  model  to  data  sets  for  two  pressure 
vessel  steels  (A515,  A533B). 

2.  Micromechanical  Constraint  Corrections 

Dodds  and  Anderson  [3,7]  show  that  by  quantifying  the  effects  of  large-scale 
yielding  (LSY)  on  the  relationship  between  microscale  crack  driving  force  (e.g.  near¬ 
tip  stresses  and  strains)  and  macroscale  crack  driving  force  (e.g.  J,  CTOD),  the  in¬ 
plane  size  effect  on  cleavage  fracture  toughness  can  be  predicted  rigorously  without 
resort  to  empirical  arguments.  Statistical  analyses  should  be  applied  only  to 
constraint  “corrected”  measures  of  fracture  toughness  [2,35].  Size  effects  become 
steadily  more  pronounced  as  load  increases  due  to  the  deviation  of  crack— tip  region 
deformations  from  the  small-scale  5delding  (SSY)  conditions  under  which  single  pa¬ 
rameter  fracture  mechanics  apply  [22,26].  When  J  (or  equivalently  CTOD)  no  longer 
describes  the  crack-tip  fields,  a  micromechanics  failure  criterion  is  required  to  es¬ 
tablish  the  near-tip  conditions  at  fracture.  Finite  element  analysis  provides  a  means 
to  quantify  the  geometry  dependent  relations  between  these  micromechanical  failure 
conditions  and  macroscale  crack  driving  force.  This  permits  (m  principle)  prediction 
of  fracture  in  any  body  from  toughness  values  measured  using  standard  specimens. 
This  section  describes  the  micromechanism  for  cleavage  fracture  developed  in  pre¬ 
vious  work  to  correct  measured  fracture  data  for  constraint  loss  due  to  large-scale 
yielding  and  then  discusses  the  basis  for  extending  the  methodology  to  include  small 
amounts  of  ductile  tearing  prior  to  cleavage  fracture. 

2. 1  Transgranular  Cleavage  Mechanism 

A  number  of  micromechanical  models  for  transgranular  cleavage  fracture  have 
been  proposed,  most  derive  from  weakest-link  statistics  (see  review  by  Wallin  [35]). 
The  weakest-link  models  assume  the  largest  or  most  favorably  oriented  fracture- 
triggering  particle  controls  the  cleavage  failure.  The  actual  trigger  event  involves  a 
local  Griffith  instability  of  a  microcrack  which  forms  at  a  microstructural  feature 
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such  as  a  carbide  or  inclusion;  satisfaction  of  of  the  Griffith  energy  balance  occurs 
when  the  critical  stress  is  reached  in  the  vicinity  of  the  microcrack  [24] .  The  size  and 
location  of  the  triggering  microstructural  feature(s)  dictate  the  fracture  toughness 
and  contribute  to  the  scatter  observed  in  results  of  cleavage  fracture  tests. 

The  Griffith  instability  criterion  implies  fracture  at  a  critical  normal  stress  near 
the  crack  tip;  the  statistical  sampling  aspect  of  the  mechanism  (i.e.,  the  probability 
of  finding  a  triggering  microfeature  near  the  crack  tip)  suggests  a  dominant  role  for 
the  volume  of  material  within  a  process-zone  over  which  the  opening  mode  stress 
exceeds  a  threshold  value  sufficient  to  initiate  cleavage.  The  probability  of  cleavage 
fracture  in  a  cracked  specimen  may  then  be  expressed  in  the  following  general  form: 

F  =  F[y(ac)]  (1) 

where  F  is  the  failure  probability  and  Vioc)  is  the  volume  of  material  ahead  of  the 
crack  tip  over  which  the  principal  stress  (ai)  exceeds  a  specified  value  Og.  This  form 
of  F  applies  to  any  fracture  process  controlled  by  maximum  principal  stress,  not  just 
weakest-link  failure  [2,35].  In  particular,  the  F  criterion  of  Eq  1  does  not  require 
material— specific  assumptions  for  the  distribution  and  strength  of  cleavage  trigger¬ 
ing  particles. 

The  present  methodology  does  not  attempt  to  predict  absolute  values  of  Jc  from 
metallurgical  parameters  that  describe  the  distribution  and  strength  of  cleavage 
triggering  particles.  Rather,  the  micromechanical  model  predicts  the  variation  of 
fracture  toughness  with  constraint  changes  for  a  given  material/temperature  by 
scaling  to  a  reference  constraint  condition.  The  crack-tip  stress  fields  in  a  test  speci¬ 
men  are  compared  to  those  for  SSY  with  T'=0.  A  J-like  parameter,  denoted  Jq,  is  ob¬ 
tained  from  this  comparison  to  the  reference  solution.  Jq  is  th®  J  to  which  the  SSY 
model  (infinite  body)  must  be  loaded  to  achieve  the  same  stressed  volume,  and  there¬ 
by  the  same  likelihood,  F,  of  cleavage  fracture,  as  in  a  finite  body.  With  this  model, 
the  assurance  of  equally  stressed  volumes  of  material  at  fracture  does  not  imply 
equivalence  of  the  individual  stress  components  beyond  SSY  conditions,  however, 
such  a  requirement  does  not  appear  strictly  essential. 


2.2  Constraint  Corrections 

Consider  the  family  of  near-tip  stress  fields  applicable  to  finite  bodies  described 
in  terms  of  the  J-Q  theory  [19,20] 
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where  Qoq  denotes  an  approximately  constant,  hydrostatic  stress  in  the  forward  sec¬ 
tor,  1 0 1  <  90°,  that  defines  the  local  deviation  of  the  finite  body  stresses  from  those 
of  the  small-scale  yielding  reference  field  (the  Q  =  0  solution  is  taken  to  be  that  for 
SSY,  T'=0).  By  employing  the  above  family  of  near-tip  stress  fields,  the  maximum 
principal  stress  also  has  the  form 
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For  any  given  value  of  Q  and  6,  oi/oq  decreases  monotonically  once  r  extends  beyond 
the  finitely  deformed  region  of  r>J/oQ.  Rearrangement  of  the  above  expression  fur¬ 
nishes  a  relation  for  the  distance  r  as  a  function  of  6  and  Oi/oq  as 
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wli6r6  §x  is  rBstrictsd  to  vflluos  ^  1  to  avoid  locations  with  tho  finitsly  doformod  ro- 
gion.  Consider  a  particular  level  of  the  principal  stress  o\laQ.  The  area  A  over  which 
the  principal  stress  exceeds  ai/  oq  is  given  by 


A  =  g  Q),  h=\^  g\{d-,  oja^,  Q)dd  . 


The  area  enclosed  by  the  contour  of  level  oiIoq  depends  on  J  as  well  as  the  triaxial- 
ity  of  the  near-tip  fields  identified  with  Q.  To  fix  ideas,  let  Aq  and  Jq  designate  the 
area  and  J  associated  with  the  Q  =  0  field,  and  let  Afb  and  JpB  designate  the  area 
and  J  associated  with  a  crack  in  a  finite  body  with  Q^O.  Then  we  have 
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Upon  initial  loading  of  the  finite  body,  Q  —  0  so  that  hps  =  Kq.  As  the  load  increases, 
plasticity  spreads  over  the  body,  Q  becomes  non-zero,  and  hpB  begins  to  deviate  from 


ho* 

For  a  given  material  and  temperature,  fracture  occurs  upon  the  attainnient  of 
equivalent  stressed  volumes  (Afb  ^  thickness  B)  for  cleavage  fracture  in  different 
specimens.  The  ratio  of  applied  J-values  in  a  finite  body  and  the  reference  Q  =  0 
stress  state  that  generate  equivalent  stressed  volumes  is  found  by  equating  areas  in 


Eqs  6  and  7  to  yield 


The  J  ratios  are  evaluated  using  Eq  8  at  each  loading  level  and  for  a  range  of 
principal  stress  values.  The  ratio  quantifies  the  size  and  geometry  dependence  of 
cleavage  fracture  toughness.  Consider,  for  example,  a  test  specimen  that  fails  at 
Jc  =  200  kJ/m^.  Suppose  the  computed  ratio  JfbUo  =  2  at  fracture  {Jfb  =  e/c)  in 
test  specimen;  then  a  much  larger  specimen  made  fi"om  the  same  material  and 
tested  at  the  same  temperature  is  predicted  to  fail  at  Jc  =  100  J/m  (the  larger  speci¬ 
men  fails  under  SSY  conditions).  Similarly,  the  fracture  toughness  ratios  for  test 
specimens  with  the  same  absolute  size  but  varying  crack  depths-to-specimen 
widths,  a/W,  may  be  quantified.  The  model  predicts  a  sharp  increase  in  fracture 
toughness  with  decreasing  alW ratio. 

An  often  overlooked  benefit  of  both  the  J-Q  description  of  near-tip  fields  and 
this  constraint  model  for  cleavage  fracture  involves  the  capability  to  model  arbitrary 
material  hardening.  A  finite  element  analysis  of  the  SSY  reference  condition  may  be 
executed  using  a  description  of  the  uniaxial  material  response  other  than  Ramberg— 
Osgood  or  the  linear,  power-law  models  (finite  strains  and  J2  flow  theory  plasticity 
may  be  modeled  as  well).  Fracture  test  specimens  (and  structural  components)  are 
then  modeled  using  the  same  procedures  with  Q  and  the  constraint  model  for  cleav¬ 
age  constructed  as  above.  These  approaches  rely  on  the  concept  of  similarity  under 
SSY  conditions;  a  concept  much  broader  than  is  inferred  by  asymptotic  analyses 
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which  invoke  incompressibility,  power-law  hardening  and  small-strain  theory  to  de¬ 
velop  expressions  for  crack— tip  fields. 

The  key  question  to  resolve  with  the  constraint  model  concerns  the  sensitivity  of 
the  Ape/Ao  and  /Jo ratios  on  o\Ioq.  Numerical  computations  [3,13,14]  reveal  a 
weak  dependence  on  OiIoq  in  common  fracture  specimens  until  the  overall  bending 
deformation  impinges  significantly  on  the  crack-tip  field.  This  occurs  more  readily 
for  low  amounts  of  strain  hardening  and  for  deeply  notched  specimens  and  much 
less  so  for  tension  dominated  geometries  and  for  materials  with  greater  amounts  of 
strain  hardening.  The  weak  dependence  of  Jpb/Jq  ratios  on  oi/oq  derives  from  the 
observation  that  stress  fields  among  fracture  specimens  vary  primarily  due  to  the 
hydrostatic  contribution  (as  described  by  the  J-Q  theory).  The  differences  in  hydro¬ 
static  stress  relative  to  the  SSY  (Q=0)  field  alter  predominantly  the  absolute  size  of 
principal  stress  contours  but  not  the  angular  form  defining  their  shapes  [8]. 

2.3  Application  of  the  Constraint  Corrections  in  Fracture  Testing 

Figure  1  provides  typical  results  obtained  through  finite  element  modeling  (plane 
strain).  The  specimens  are  SE(B)s  with  atW  =  0.1,  0.5  modeled  with  a  Ramberg-Os- 
good  material  having  strain  hardening  exponents  n  =  5, 10.  Values  of  JpB  and  Jo  are 
plotted  on  separate  axes  to  facilitate  the  constraint  correction  of  experimental  data. 
Points  on  the  curves  describe  (Jfb,  Jq)  pairs  that  produce  equal  stressed  volumes  of 
material  in  the  test  specimen  and  in  SSY.  Path  A-B-C  indicated  in  Fig.  1  (b)  illus¬ 
trates  the  procedure  to  correct  for  the  effects  of  large-scale  yielding  in  measured 
cleavage  fracture  toughness  data  (Jc  value  at  A)  by  determining  the  geometry  inde¬ 
pendent  fracture  toughness  (Jq  value  at  C).  To  simplify  application  of  the  constraint 
correction  procedures,  curve  fits  to  finite  element  results  of  the  form  shown  in  Fig.  1 
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FIG.  1 — Cleavage  constraint  models  for  stationary  cracks  in  SE(B)  specimens 
with  a/W=0.1,  0.5;  n  =  5,  10;  E/oq  =  500;  using  contour  areas  for  oiloo=  2.7. 
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are  available  for  a  wide  range  of  hardening  exponents  and  alW  ratios  [8].  Moreover, 
crack  opening  stresses  on  the  remaining  ligament  may  be  used  to  compute  the 
JfbUq  ratios  rather  than  the  more  complex  contour  area  procedure  [8].  Applications 
of  this  technique  to  structural  steels,  including  A36,  A515  and  A533B  [12,29,35], 
have  been  successful  in  removing  the  geometry  dependence  of  Jc  values. 

Three-dimensional  analyses  in-progress  for  the  same  SE(B)  specimens  (with 
BxB  cross  sections)  reveal  significant  differences  with  the  plane-strain  correction 
model  for  the  deep-notch  specimen  but  essentially  no  differences  with  plane-strain 
results  for  the  shallow-notch  specimen.  The  3-D  analyses  show  that  SSY  conditions 
for  the  a/W  =  0.5  specimen  are  maintained  to  greater  deformation  levels  (aao/J  =  50 
vs.  150  for  plane-strain)  on  the  mid-thickness  plane;  at  positions  on  the  crack  front 
mid— way  between  the  center— plane  and  outside  surface,  the  constraint  levels  fall  be¬ 
low  those  predicted  by  the  plane-strain  model.  For  the  shallow-notch  specimen, 
constraint  conditions  on  the  mid-thickness  plane  closely  match  those  predicted  by 
the  plane-strain  model.  Introduction  of  the  finite  thickness,  B,  requires  modification 
of  the  constraint  correction  model  to  accommodate  the  interaction  of  thickness  with 
a/W  ratio,  W/B  ratio  and  material  strain  hardening.  The  concept  of  an  effective 
thickness,  Bgff,  is  being  explored  to  describe  constraint  across  the  crack  front  relative 
to  the  mid— thickness  constraint  in  a  meaningful  way  for  cleavage  fracture.  The  rela¬ 
tive  effects  of  crack  growth  in  a  plane-strain  vs.  a  fully  3-D  model  (with  curved 
crack  fronts  during  growth  allowed)  appear  to  be  a  very  important,  yet  currently  un¬ 
resolved  issue  (the  computational  resources  and  modeling  complexity  for  growth 
analyses  in  3-D  are  much  greater). 

3.4  Extension  to  Include  Effects  of  Ductile  Tearing 

Small  amounts  of  ductile  crack  extension,  typically  less  than  1-2  mm,  preceding 
cleavage  fracture  may  alter  the  structure  and  severity  of  the  near-tip  stress  fields. 
The  extending  tip  develops  a  sharp  opening  profile,  relative  to  the  blunt  profile  at 
initiation  of  tearing,  which  elevates  the  stress  and  decreases  the  distance  from  the 
crack  tip  to  the  peak  stress  location.  When  SSY  conditions  prevail,  stresses  at  dis¬ 
tances  greater  than  J /oq  from  the  tip  increase  as  well  due  to  additional  elastic 
strains  as  the  plastic  zone  pushes  forward  into  unyielded  material.  Crack  growth 
also  subjects  a  larger  volume  of  material,  AY  « (J /oq)  x  Aa  x  B,  to  high  stress  levels 
near  the  tip  and  thus  raises  the  probability  of  sampling  a  critical,  firacture— trigger¬ 
ing  particle. 

At  the  metallurgical  length-scale  characteristic  of  carbides  and  inclusions,  tear¬ 
ing  and  large  plastic  deformation  may  affect  the  fundamental  cohesive  strength 
across  cleavage  planes.  More  likely  at  this  length— scale,  tearing  makes  the  crack 
front  considerably  more  torturous  than  the  blunted  tip  of  the  stationary  crack,  thus 
increasing  the  cleavage  energy  needed  to  initiate  a  macroscopic  fracture.  In  addition, 
material  damage  caused  by  the  formation  and  growth  of  voids  within  the  zone  of  fi¬ 
nite  strains  may  alter  somewhat  the  macroscopic  stress  field  [15].  These  are  unre¬ 
solved  issues. 

To  extend  the  constraint  correction  model  given  present  knowledge,  these  metal¬ 
lurgical  aspects  of  tearing  are  considered  to  have  secondary  importance  relative  to 
the  increased  near— tip  stresses  and  the  additional  volume  of  material  subjected  to 
higher  stresses  as  a  consequence  of  the  crack  extension.  The  effects  of  ductile  tearing 
thus  appear  as  changes  in  the  stressed  volumes  of  material,  XjB,  where  crack 
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growth  changes  the  area  enclosed  within  principal  stress  contours  relative  to  a  sta¬ 
tionary  crack  loaded  at  the  same  JpB. 

Figure  2  illustrates  schematically  the  impact  of  tearing  on  the  constraint  correc¬ 
tion  model.  The  specimen  follows  tte  model  for  a  stationary  crack  until  the  onset  of 
ductile  tearing  at  JpB  =  Jic-  Under  continued  loading,  the  new  model  including 
growth  predicts  an  increase  of  «/o-values  relative  to  continued  loading  of  the  station¬ 
ary  crack,  i.e.,  tearing  increases  the  constraint.  The  increase  of  Jo-values  due  to 
tearing  depends  on  the  amount  of  plastic  deformation  in  the  specimen  at  J/c  and  on 
the  iJ-curve  followed  prior  to  cleavage  fracture.  To  simplify  initial  development  of 
the  model,  a  constant  tearing  modulus,  Tj  =  (E/ol)dJ /da,  is  adopted  to  describe 
the  R-curve  for  t5^ical  extensions  of  1-2  mm.  Consequently,  the  evolution  of  crack- 
tip  constraint  in  simple  fracture  specimens  with  ductile  growth  may  be  expressed  by 
a  function  of  the  form 

^FB  _  ip/^i  „  E  ^0  ^'^0  Aa  rp\  (9) 

Jo  W’  Jj^  ’ 

where  C  denotes  the  “constraint”  function  for  cleavage  fracture  with  prior  ductile 
tearing.  The  non-dimensional  groups  Woq/Jj^,  aJW,  Aa/W  and  Tj  describe  the  de¬ 
pendence  on  crack  growth  and  absolute  size,  while  E/oq  and  n  define  the  material 
5deld  and  hardening  properties.  Although  included  Eq  9,  a  weak  dependence  of  C  on 
OiIoq  is  expected. 

Jq— values  derived  from  both  instantaneous  and  cumulative  contour  areas  are 
provided  for  initial  development  of  the  growth  correction  model  (see  Fig.  3).  The 
instantaneous  areas  enclosed  within  principal  stress  contours  are  computed  from  a 
snapshot  of  the  stress  field  ahead  of  the  growing  crack —  points  on  such  contours  all 
lie  within  the  forward  sector,  1 6 1  <  90°  (see  area  X  in  Fig.  3).  Figure  3  also  illus¬ 
trates  the  construction  of  cumulative  contour  areas  using  the  envelope  of  maximum 
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principal  stresses  experienced  by  material  along  the  crack  plane  during  growth.  The 
cumulative  areas,  defined  as  areas  enclosed  within  the  contour  envelope,  may  be 
substantially  larger  than  instantaneous  areas  which  leads  to  the  larger  Jo-values 
indicated  in  Fig.  2  for  the  cumulative  area  correction. 

The  proposed  constraint  model  including  crack  growth  maintains  the  relative 
simpKcity  of  the  original  version  which  corrects  only  for  (in-plane)  effects  of  large- 
scale  yielding.  Furthermore,  development  of  the  new  model  to  account  only  for  tear¬ 
ing  effects  on  the  microscale  driving  force  (the  local  crack  tip  stresses  and  stressed 
volume),  implies  that  corrected  toughness  values,  with  and  without  ductile  tearing 
prior  to  cleavage,  should  receive  equal  weight  in  statistical  analyses,  e.g.,  Weibull 
models  [34,  35]. 

3.  Numerical  Models  for  Growth  Analyses 

Finite  element  models  of  the  type  shown  in  Fig.  4  are  employed  to  perform  plane 
strain,  nonlinear  analyses  of  SE(B)  specimens  including  small  amounts  of  ductile 
crack  growth.  Both  shallow  crack  (a/W=0.1)  and  deep  crack  (a/W=0.5)  specimens  are 
analyzed.  Symmetry  about  the  crack  plane  permits  modeling  of  only  one-half  of  each 
SE(B)  specimen.  The  finite  elements  are  conventional  eight-noded  isoparametrics 
with  reduced  (2x2)  Gauss  integration.  Square  elements  in  the  crack-tip  region  and 
along  the  crack  plane  are  defined  to  simplify  the  growth  analyses  and  to  permit  uni¬ 
form  increments  of  crack  extension.  Crack  opening  profiles  and  stresses  ahead  of  the 
growing  crack  seem  to  be  resolved  adequately  when  the  crack— tip  element  size,  Lg,  is 
smaller  than  «0.5  X  the  CTOD  at  initiation  (djc).  Stationary  crack  solutions  for  com¬ 
parison  are  also  generated  with  these  models. 

The  finite  element  models  are  loaded  by  displacement  increments  imposed  on  a 
single  centerplane  node  as  shown  in  Fig.  4.  An  improved  form  of  the  traditional  node 
release  technique  is  employed  to  implement  crack  extension.  The  imposed  load- 
point  displacement  is  adjusted  iteratively,  concurrent  with  a  one-element  crack  ex- 


FIG.  Z— Schematic  illustrating  construction  of  principal  stress  envelope  for 
growing  cracks. 
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Finite  element  mesh  SE(B)  specimen  with  deep  crack 
(a/W  =  0.5).  Elements  are  8-node  isoparametrics. 


Le/W  =  0.001 


Crack- plane  elements  are  square  to  simplify  crack 
(b)  growth  analyses.  Lq  denotes  side  length  of  the  square 
crack  plane  elements. 

FIG.  4 — Finite  element  models  for  SE(B)  specimens. 


tension,  such  that  the  crack  growth  criterion  is  always  satisfied  at  completion  of  the 
node  release.  This  procedure  retains  the  computational  convenience  of  the  node  re¬ 
lease  technique  while  eliminating  convergence  problems  due  to  artificial  unloadings 
encountered  with  the  traditional  implementation  (see  [9]  for  additional  details). 

The  uniaxial  stress-strain  relationship  follows  the  Ramberg-Osgood  model, 
s/sq  =  o/oq  -F  a(a/ao)”-  Two  t5q)ical  values  of  the  material  hardening  coefficient 
(n  =  5  and  10)  are  selected  in  the  analyses  of  SE(B)  specimens  to  represent  high 
(n  =  5)  and  medium(n  =  10)  hardening  materials.  All  computations  use  E/ao  =  500, 
Poisson’s  ratio  v  =  0.3  and  a  =  1. 

The  computations  reported  here  are  carried  out  using  a  large-rotation,  finite 
strain  constitutive  model  based  on  the  incremental,  flow  theory  of  plasticity  [9,10\.  A 
companion  study  using  a  conventional,  small-strain  (incremental)  plasticity  for¬ 
mulation  predicts  essentially  identical  results  for  crack  opening  profiles  and  stresses 
ahead  of  the  growing  tip.  Recent  work  by  Liu  and  Drugan  reach  similar  conclusions 
[23]. 
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3. 1  Crack  Growth  Parameters 

The  specified  value  of  Jjc,  relative  to  specimen  geometry,  size  and  material  flow 
properties,  sets  the  severity  of  plastic  deformation  in  an  SE(B)  specimen  at  the  initi¬ 
ation  of  ductile  growth.  Values  of  Who/J/c  =  480^80  are  defined  to  achieve  conditions 
at  initiation  of  growth  ranging  from  essentially  SSY  to  fully  formed  plastic  hinges. 

Once  J  reaches  the  specified  initiation  toughness,  Jjc,  subsequent  crack  exten¬ 
sion  follows  an  jR-curve  defined  by  a  constant  tearing  modulus,  Tj.  Selected  values 
for  the  tearing  modulus  of  15,  40  and  83  represent  relatively  low  values  characteris¬ 
tic  of  higher-strength  steels  with  moderate-to-low  strain  hardening.  Values  on  the 
order  of  200-400  occur  during  the  early  stages  of  crack  growth  in  some  lower- 
strength  steels  with  moderate-to-high  strain  hardening.  These  large  values  ap¬ 
proach  one-half  the  slope  of  the  blunting  line,  dJ/da  =  2ao  (defines  a  tearing 
modulus  of  Tjbl  =  /oq).  The  lower  values  of  Tj  specified  for  the  SE(B)  analyses 

accentuate  the  effects  of  crack  growth  —  much  larger  values  of  tearing  modulus 
drive  the  response  toward  the  stationary  crack  solution,  i.e.,  the  loss  of  constraint 
due  to  large  increases  in  J  and  plastic  deformation  overwhelm  the  increased 
constraint  from  small  crack  extensions.  Analyses  for  materials  with  high  tearing 
modulus  (200-500)  are  also  performed  to  support  application  of  the  growth 
constraint  model  to  currently  available  experimental  results. 

3.2  J-Values  During  Crack  Growth 

J-values  during  crack  growth  analyses  in  the  SE(B)  specimens  are  computed 
with  the  domain-integral  procedure  [18].  Figure  5  shows  the  typical  variation  of  J- 
values  with  increasing  generalized  radius  of  the  integration  domain  for  both  a  sta¬ 
tionary  crack  and  a  growing  crack.  A  domain  “ring”  denotes  a  sequence  of  connected 
elements,  enclosing  the  tip  at  a  remote  distance,  which  make  non-zero  contributions 
to  J.  Starting  with  the  ring  of  elements  at  the  current  crack  tip,  J  is  computed  inde¬ 
pendently  over  each  successive  ring  of  elements.  Figure  5  (a)  shows  a  weak  path  de¬ 
pendence  of  J  on  the  domain  size  (i.e.,  distance  of  ring  from  tip)  for  the  stationary 
crack  loaded  to  J  =  J/c-  The  small  degree  of  path  dependence  just  outside  the  crack- 
tip  elements  arises  from  nonproportional  straining  in  the  incremental  plasticity 
theory.  Once  the  domains  extend  beyond  10%  of  the  remaining  ligament,  J  becomes 
path-independent;  this  value  is  referred  to  here  as  the  far-field  J  and  is  denoted 

Jfar- 

In  contrast  to  the  stationary  crack,  J-values  for  the  growing  crack  exhibit  a 
much  stronger  path— dependence  as  shown  in  Fig.  5  (b).  Near  zero  J— values  are  com¬ 
puted  over  domains  very  close  to  the  crack— tip.  Figure  5  (b)  shows  that  the  J— value 
at  Aa  =  10  X  djc  reaches  a  nearly  path-independent  value  once  domains  extend  be¬ 
yond  20-25%  of  the  remaining  ligament.  The  necessity  for  J  evaluation  at  large  dis¬ 
tances  from  the  growing  crack  limits  the  maximum  amount  of  crack  growth  per¬ 
mitted  in  the  SE(B)  analyses.  As  the  largest  domain  is  redefined  for  each  increment 
of  growth,  the  domain  may  extend  into  the  more  coarse  mesh  remote  from  the  tip  or 
reach  the  specimen  edges  ahead  of  or  behind  the  crack  tip. 

In  the  experimental  determination  of  a  J/j— Aa  response,  J-values  are  estimated 
as  a  fraction  of  elastic  and  plastic  work  done  on  the  specimen.  By  treating  the  finite 
element  P-Ar  j.n  results  (including  crack  growth)  as  experimental  records,  a  compar¬ 
ison  of  the  domain  integral  and  experimental  J  definitions  becomes  possible.  For 
deep  notch  SE(B)  specimens,  ASTM  E 1152-87  [4]  outlines  the  procedure  to  estimate 
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JUlC  J/Jlc 


Domain  Radius  /  (W-a)  Domain  Radius  /  (IV-ao) 

(a)  Stationary  Crack  (b)  Growing  Crack 

FIG.  5 — J-integral  variation  over  domain  rings  for  both  stationary  and 
growing  cracks  for  SE(B)  with  a/W=0.5. 

J-values  (Jastm)  from  experimentally  measured  load,  load-line  displacement  and 
crack  length  records.  (The  ASTM  procedure  assumes  nonlinear  elastic  material  re¬ 
sponse  as  well  and  thus  the  Jastm  values  are  often  referred  to  as  “deformation”  J- 
values,  Jd).  For  both  hardening  exponents  (n  =  5  and  10),  the  finite  element  Jfar  and 
Jastm  values  are  indistinguishable  over  the  full  crack  growth  history.  Such  close 
agreement  of  finite  element  and  experimental  estimates  for  J-values  with  crack 
growth  is  essential  for  subsequent  application  of  the  growth  constraint  model. 

4.  Growth  Effects  on  Crack-Tip  Fields 

4. 1  Crack  Opening  Profiles 

Figure  6  compares  crack  opening  profiles  for  stationary  crack  and  growing  cracks 
in  SE(B)  specimens.  The  effects  of  different  tearing  modulus  are  shown  for  crack 
growth  initiated  at  the  same  normalized  value  of  J/c .  Profiles  for  growing  cracks  ex¬ 
hibit  a  macroscopically  “sharp”  tip  for  all  values  of  tearing  modulus,  in  contrast  to 
the  severely  blunted  tip  of  the  stationary  crack.  While  the  different  values  of  tearing 
modulus  strongly  affect  the  immediate  crack  opening  angle  (CTOA),  as  indicated  on 
the  figure,  there  is  negligible  effect  upon  further  opening  behind  the  initial  (blunted) 
tip  when  displacements  are  normalized  by  J Ioq.  Following  a  short  transient  period 
when  the  crack  grows  through  the  first  element  at  the  stationary  tip,  the  crack-tip 
opening  angles  remain  constant  and  increase  with  the  specified  Tj.  The  constancy  of 
CTOA  with  Tj  observed  previously  in  SSY  analyses  of  crack  growth  [25,28]  carries 
forward  to  the  SE(B)  specimens. 

4.2  Crack  Opening  Stresses  on  Remaining  Ligament 

Figure  7  illustrates  the  typical  effects  of  growth  on  the  crack  opening  stress 
which  acts  on  the  remaining  ligament.  The  analyses  consider  two  geometries 
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(a/ W  =  0.1,  0.5)  and  two  hardening  exponents  {n  =  5, 10),  with  results  for  each  com¬ 
bination  of  these  parameters  shown  for  one  value  of  the  tearing  modulus  {Tj  =  83). 
Crack  growth  begins  in  the  analyses  at  the  common  value  Jfc  =  Wao/240.  The  results 
for  each  case  contain  two  curves:  one  for  the  growing  crack  of  length  ao  +  Aa  and  one 
for  the  stationary  crack  of  length  ao,  with  both  configurations  loaded  to  the  same  J. 
The  amount  of  crack  extension  relative  to  the  CTOD  at  initiation  varies  with  the 
strain  hardening  even  though  the  final  J-values  are  1.8e7/c  in  each  case.  For  refer¬ 
ence,  each  figure  also  includes  the  corresponding  steady-  growth  solution  in  SSY 
given  by  Varias  and  Shih  [31]  . 

Sharpening  of  the  crack  tip  during  growth  shifts  the  location  of  the  peak  stress 
nearer  to  the  tip,  with  the  largest  (relative)  shift  for  the  lower  amount  of  strain 
hardening.  The  peak  stress  location  for  both  high  and  moderate  hardening  converges 
to  r«0.25  J/ao.  As  expected,  the  shallow  crack  specimens  reveal  significant  loss  of 
constraint  relative  to  the  deep  crack  specimens,  with  stresses  for  n  =  10  showing  the 
most  differences  [compare  stresses  at  the  same  rl(JI  (Tq)]  •  At  distances  from  the  tip 
greater  than  about  J/ao,  crack  growth  in  the  SE(B)  specimens  does  not  alter  signifi¬ 
cantly  stresses  on  the  remaining  ligament.  This  behavior  contrasts  sharply  with 
stresses  for  the  steadily-growing  crack  in  SSY;  stresses  larger  than  those  for  the 
non— growing  crack  are  predicted  over  much  greater  distances  from  the  tip  [31].  In 
these  SE(B)  specimens,  crack  growth  initiates  and  grows  under  conditions  of  large- 
scale  5delding  on  the  ligament  which  appears  to  suppress  other  mechanisms  of  stress 
elevation  beyond  the  sharpening  effect. 

For  the  n  =  10  material,  the  amplitude  of  peak  stress  during  crack  growth  ex¬ 
ceeds  the  amplitude  of  the  corresponding  stationary  crack  stress  and,  for  the  deep 
notch  specimen,  the  peak  stress  reaches  the  steady  growth  value  in  SSY.  For  the 
high  hardening  material  (n=5),  the  growth  stresses  appear  to  fall  below  the  peak 
values  without  growth.  The  relative  coarseness  of  the  element  mesh  is  clearly  evi¬ 
dent  in  the  stress  fields  without  growth  for  the  n  =  5  material;  the  crack— tip  opening 
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FIG.  6 — Opening  profiles  for  stationary  and  growing  cracks  for  shallow  and  deep 
notch  SE(B)  specimens. 
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displacement  for  this  case  is  comparable  to  the  tip  element  size.  We  expect  that  re¬ 
analysis  using  a  finer  mesh  at  the  tip  for  n  =  5  would  predict  somewhat  lower  peak 
stresses  for  the  stationary  crack.  Comparisons  of  crack  opening  stresses  computed 
for  other  cases  reveal  that  variations  of  Jjc  and  the  strain  hardening  exponent  affect 
the  stresses  more  significantly  than  variations  of  the  tearing  modulus  when  dis¬ 
tances  are  scaled  by  JIoq  as  in  Fig.  7. 

4.3  Principal  Stress  Contours 

Crack  opening  stresses  on  the  remaining  ligament  may  not  reflect  changes  in 
stresses  away  from  the  crack  plane  due  to  growth.  Figure  8  (a-b)  show  the  effects  of 
crack  growth  on  the  instantaneous  principal  stress  contours  for  shallow  and  deep 
crack  SE(B)  specimens.  Crack  growth  begins  at  the  common  value  Jjc  =  WoqI2A0 
with  results  given  in  the  figure  for  Tj  =  83  and  n  =  10.  Comparisons  for  contours  oth¬ 
er  than  cTi/ao  =  2.7  shown  in  the  figure  reveal  similar  trends.  In  each  case  the  sta¬ 
tionary  crack  and  growing  crack  contours  correspond  to  the  same  applied  J. 

By  using  the  current  value  of  J/oq  to  normalize  distances  from  the  crack  tip  on 
plots  of  principal  stress  contours,  deviations  in  the  stress  triaxiality  due  to  crack 
growth  and  increased  loading  (larger  J  value)  become  very  clear.  Under  SSY  condi¬ 
tions  the  contours  for  a  specified  oi/oq  ratio  coincide  at  all  J-values  when  expressed 
on  these  normalized  axes.  When  contours  in  a  specimen  at  different  loads  coincide, 
the  near-tip  constraint  remains  unchanged  (the  absolute  size  of  the  contour  in¬ 
creases  with  J  but  at  the  same  rate  for  the  two  loading  levels).  When  constraint  in¬ 
creases  (decreases)  under  increased  loading  or  crack  growth,  the  normalized  con¬ 
tours  must  also  increase  (decrease)  in  size. 

The  normalized  contours  for  the  growing  crack  are  distinctly  larger  than  the  con¬ 
tours  for  the  stationary  crack,  especially  for  the  shallow  crack  SE(B)  specimen.  For 
the  a/W  =  0.5  specimen,  the  opening  mode  stresses  on  the  remaining  ligament  do 
not  indicate  correctly  the  effects  of  crack  growth  on  stresses  elsewhere  in  the  for¬ 
ward  sector.  The  contours  shown  in  the  figure  are  those  at  the  maximum  crack 
growth  considered  of  « 10  x  contours  examined  at  less  crack  growth  show  some¬ 
what  larger  differences. 

Figure  8  (c-d)  summarize  the  areas  enclosed  by  the  instantaneous  contours  for  a 
range  of  0\Ioq  ratios.  Areas  are  normalized  by  the  similarity  factor  (in  SSY 
with  Aa  =  0,  the  quantity  remains  constant  under  increased  loading). 

Three  sets  of  results  for  comparison  are  shown  on  each  figure:  (1)  areas  for  SSY 
without  crack  growth,  (2)  areas  for  the  SE(B)  specimen  at  the  maximum  amount  of 
crack  extension  considered,  and  (3)  areas  for  the  SE(B)  specimen  without  crack 
growth  hut  loaded  to  the  same  J-value  as  the  growing  crack.  Figure  8  (c)  for  the 
shallow  crack  demonstrates  the  dramatic  loss  of  constraint  caused  by  the  extensive 
plastic  deformation;  the  areas  differ  by  a  factor  of  =:  100  from  SSY  for  o\loi  =  2.7. 
Crack  growth  clearly  increases  the  constraint,  elevating  the  areas  by  factors  of  2-3 
above  stationary  crack  levels.  The  deep  crack  specimens  maintain  much  higher 
constreiint  as  evidenced  by  the  smaller  differences  in  contour  areas  compared  to  SSY. 
Crack  growth  elevates  near-tip  stress  triaxiality  and  increases  areas  enclosed  with¬ 
in  the  principal  stress  contours.  However,  the  effects  of  growth  on  the  contours  are 
much  less  pronounced  than  for  the  shallow  crack. 

Not  surprisingly,  the  differences  in  contour  areas  between  the  stationary  and 
growing  crack  increases  with  the  OiIoq  ratio,  that  is,  contours  near  the  growing  tip 
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FIG.  1— Effects  of  growth  on  crack  opening  stress  in  SE(B)  specimens.  For  each 
a/W,  growing  cracks  and  stationary  cracks  are  loaded  to  the  same  J. 
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SE(B)  specimens.  For  each  alW,  growing  cracks  and  stationary  cracks  are  loaded  to 


the  same  J. 
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sense  the  increased  peak  stress  caused  by  the  sharpening.  The  trends  in  contour 
areas  shown  in  these  figures  are  entirely  consistent  with  observations  about  growth 
effects  on  constraint  made  previously  on  the  basis  of  crack  opening  stresses  acting 
on  the  remaining  ligament. 

In  summary,  our  extensive  studies  [30]  of  SE(B)  specimens  with  crack  growth  to 
examine  the  effects  of  n,  Jjc,  Tj  reveal  a  complex  interaction  of  constraint  trade-offs. 
Increased  strain  hardening  increases  the  level  of  constraint  and  reduces  the  effect  of 
growth  on  the  near  tip  fields.  Larger  values  of  J/c  increase  the  impact  of  crack 
growth  on  the  near-tip  stresses — crack  growth  has  more  effect  when  the  stress 
fields  at  initiation  have  lost  significant  constraint  relative  to  SSY.  Larger  values  of 
the  tearing  modulus  reduce  the  net  change  in  the  constraint  due  to  growth — a  large 
teeiring  modulus  requires  large  J  increases  to  extend  the  crack  which  in-turn  accel¬ 
erate  the  constraint  loss  by  additional  large-scale  yielding.  Consequently,  crack 
growth  has  more  effect  as  7>->0  and  negligible  effect  as  Tj->Tjbl  (the  blunting  line 
slope). 

5.  Tearing  Effects  on  Constraint  Model 

The  extensive  finite  element  analyses  of  SE(B)  specimens  that  include  the  effects 
of  crack  growth  provide  a  basis  to  extend  the  previously  developed  constraint  model 
for  cleavage  fracture.  The  SSY  stress  fields  (without  crack  growth)  again  define  a 
convenient  reference  condition  to  assess  constraint  effects.  The  JpbUq  ratios  are 
computed  using  Eq  8  for  two  afW  ratios  (0.1,  0.5),  for  two  strain  hardening  expo¬ 
nents  (n=5,  10)  and  for  a  wide  range  of  values  for  J/c  and  tearing  modulus.  The  com¬ 
putations  provide  JfbUq  ratios  based  on  areas  within  the  instantaneous  principal 
stress  contours  and  on  cumulative  areas  within  the  envelope  of  contours  constructed 
over  the  history  of  crack  growth. 

Figures  9  and  10  show  the  typical  variation  of  JfbUq  ratios  with  principal  stress 
contour  found  in  the  parameter  study.  These  figures  illustrate  the  effects  of  alW  ra¬ 
tio  and  hardening  exponent  for  a  common  value  of  J/c  and  Tj,  and  also  compare  the 
JfbUo  ratios  for  instantaneous  and  cumulative  areas.  The  solid  lines  denote  JfbUq 
ratios  at  the  onset  of  crack  growth,  JpB  =  Jlc  \  each  dotted  line  indicates  the  ratios  at 
a  fixed  amount  of  crack  growth.  When  the  JfbUq  ratio  increases,  constraint  (stress 
triaxiality)  in  the  specimen  decreases. 

Consider  first  the  shallow  crack  results  shown  in  Fig.  9.  The  significant  features 
include:  (1)  high  strain  hardening  (n  =  5)  limits  constraint  loss;  (2)  JfbUq  ratios  for 
crack  growth  are  less  dependent  on  the  principal  stress  ratio  than  are  the  ratios  for 
stationary  cracks  (for  n=10,  the  stationary  crack  has  lost  similarity  of  the  principal 
stress  contours  at  e7/c);  (3)  JpB^^O  ratios  based  on  cumulative  areas  are  smaller  than 
ratios  based  on  instantaneous  areas,  and  (4)  most  importantly,  JfbUq  ratios  de¬ 
crease  relative  to  the  stationary  crack  values  with  increased  crack  growth  and  in¬ 
creased  loading.  This  last  feature  contrasts  sharply  with  the  behavior  of  stationary 
cracks  which  show  continually  increasing  JfbI^o  ratios  under  additional  loading. 
Constraint  in  the  shallow  crack  specimen  with  n  =  10  increases  continuously  with 
crack  growth.  The  response  in  Fig.  9  (a)  for  n  =  5  computed  using  instantaneous 
areas  provides  a  particularly  interesting  result  in  that  the  constraint  remains  essen¬ 
tially  constant — the  increased  constraint  caused  by  crack  growth  cancels  almost  ex¬ 
actly  the  constraint  loss  due  to  additional  plastic  deformation  under  continued  load¬ 
ing. 
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FIG.  9 — Effects  of  principal  stress  ratio  and  instantaneous  vs.  cumulative  areas 
on  Jo  for  a  shallow-notched  SE(B)  specimen  (alYf-0.1). 
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Consider  next  the  deep  crack  results  shown  in  Fig.  10.  The  overall  trends  remain 
similar  except:  (1)  JpbUq  ratios  at  the  onset  of  tearing  are  weakly  dependent  on  the 
principal  stress  ratio  for  both  low  and  high  amounts  of  hardening  and  (2)  J/rjg/Jo  ra¬ 
tios  for  growth  based  on  cumulative  area  for  n  =  10  show  a  more  noticeable  depen¬ 
dence  on  the  principal  stress  contour.  During  crack  growth  with  the  corresponding 
increase  in  J,  the  specimen  deformation  increases  substantially  and  the  bending 
field  impinges  on  the  crack-tip  region. 

Figure  11  recasts  these  Jpb/Jq  ratios  for  oi/oq  =  2.7  into  the  usual  format  of  the 
constraint  model.  Results  of  similar  computations  for  additional  values  of  J/c  are  giv¬ 
en  as  well.  Solid  lines  denote  the  stationary  crack  response  while  dashed  lines  de¬ 
note  the  growth  response  in  terms  of  the  instantaneous  and  cumulative  contour 
areas.  Each  curve  for  the  specific  a/W  ratio  terminates  at  the  same  amount  of  crack 
growth  in  the  specimen  iAomax  is  indicated  for  each  set  of  results).  A  reference  line 
indicating  AJq  =  AJpB  is  provided  on  the  figures  to  aid  in  assessing  the  relative 
constraint  changes. 


On  these  and  subsequent  figures,  we  normalize  J-quantities  using  the  flow 
stress  (ofiow)  rather  than  oq  to  make  the  results  less  sensitive  to  the  actual  estimate 
of  the  hardening  exponent.  For  a  Ramberg-Osgood  material  representation  used  in 
the  finite  element  computations,  an  estimate  for  the  flow  stress  is  given  by 


(iVAo)« 

exp(iV) 


(10) 


where  N  =  1/n.  The  remaining  ligament  prior  to  crack  growth,  b  =  W-ao,  provides  a 
suitable  specimen  dimension  for  normalizing  J  in  the  constraint  model. 


This  figure  shows  clearly  the  effects  of  crack  growth  on  the  effective  J  for  cleav¬ 
age  fracture  (Jq)  for  increasing  Jjc  at  a  common  value  of  the  tearing  modulus 
(Tj=40).  All  dashed  lines  in  Fig.  11  lie  above  the  solid  lines,  which  indicate 
constraint  increases  with  crack  growth.  Significant  observations  from  this  figure  in¬ 
clude:  (1)  the  constraint  model  using  the  cumulative  areas  always  predicts  larger 
values  of  effective  fracture  driving  force  (Jq);  (2)  the  greater  impact  of  growth  for  the 
shallow  crack  specimens  is  apparent  by  the  much  sharper  deviation  of  the  growth 
curves  from  the  stationary  crack  curve;  (3)  for  this  combination  of  material  parame¬ 
ters,  the  cumulative  growth  curves  have  slopes  of  =1,  i.e.,  increments  of  JpB  applied 
to  the  specimen  after  the  onset  of  growth  are  transmitted  fully  into  increments  of 
increased  driving  force  for  cleavage. 


5. 1  Normalizing  Jjc  Effects 

Figure  11,  and  others  constructed  for  a  range  of  tearing  modulus  values,  provides 
the  motivation  for  development  of  a  normalization  procedure  to  minimize  the  depen¬ 
dence  of  the  model  on  values  of  Jjc-  The  family  of  growth  constraint  curves  for  differ¬ 
ent  values  of  the  tearing  modulus  (at  a  common  Jjc)  simply  translates  with  in¬ 
creased  values  for  Jjc  while  retaining  a  nearly  constant  orientation  relative  to  the 
gradient  of  the  constraint  curve  for  the  stationary  crack.  Figure  12  (a-b)  demon¬ 
strate  the  effects  of  J/c  on  the  growth  portions  of  the  constraint  model  when  viewed 
in  this  framework.  In  Fig.  12  (a),  m  represents  the  gradient  of  the  stationary-crack 
constraint  model  a.tJic,m-  dJQ/dJpB;  a  denotes  the  angle  change  between  the  tan¬ 
gent  line  to  the  stationary  curve  (which  has  slope  m)  and  a  similar  line  drawn  tan¬ 
gent  to  the  growth  curve  at  the  same  J/c  value.  When  the  J/c  value  increases  from 
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c/q/ i^^low)  X  1  0^  w/q/ i^^low)  X  1  0^ 


UpBlibOfiow)  ^FB/ibOfiow) 

(a)  a/W  =  0.1  (b)  a/W  =  0.5 


FIG.  11 — Comparison  of  constraint  corrections  for  growth  based  on  areas 
enclosed  within  instantaneous  and  cumulative  principal  stress  contours  for  SE(B) 
specimens  (Tj  =  40,  oiIoq-  2. 7). 


to  m  decreases  significantly  from  m\  to  m2  as  indicated  on  the  figure.  Figure 
13*^ provides  values  of  m  for  SE(B)  specimens  containing  shallow  and  deep  cracks  for 
hardening  exponents  n  =  5, 10. 

Figure  12  (b)  shows  growth  curves  defined  at  increasing  values  of  Jjc  after  the 
applicable  Jjc  value  is  subtracted  from  the  total  values  of  JpB  and  Jq.  The  shifted 
curves  now  show  only  the  increments  of  J  that  occur  during  growth,  denoted 
and  AJo.  and  defined  by: 


~  ^FB  ^Ic 
AJq  =  Jq  —  J^Q 


(11) 

(12) 


where,  for  the  remainder  of  this  section,  all  J-quantities  are  implicitly  normalized 
by  the  applicable  bafiow'y  Jq  represents  the  value  of  Jo  at  JpB  =  die  -  Using  these  new 
quantities,  the  nearly  equal  values  of  a  at  the  different  values  of  J/c  become  more 
apparent.  The  a  values  exhibit  a  weak  dependence  on  J/c  for  the  alW=  0.5  specimen, 
but  are  strongly  independent  of  J/c  for  the  a/W  =  0.1  specimen.  Since  a  remains  rel¬ 
atively  constant,  the  dependence  of  the  growth  curves  with  J/c  may  be  approximated 
as  a  linear  function  of  the  gradients  (m)  for  the  stationary  crack. 


A  simple  procedure  now  defines  the  relationship  between  Jo  values  for  two  differ¬ 
ent  J/c  values.  Consider  points  A  and  B  in  Fig.  12  (a)  which  represent  specimens 
with  the  same  Tj  and  same  (normalized)  amount  of  crack  growth  but  with  different 
(normalized)  J/c  values  —  configurations  (1)  and  (2)  thus  may  represent  specimens 
differing  in  absolute  size  but  having  the  same  absolute  value  of  J/c.  The  constraint 
corrected  values  for  these  two  points  are  labelled  Jq  and  J®  in  Fig.  12  (a).  During 
crack  growth  along  (2)  to  reach  point  B,  JpB  changes  by  A  J^^  =  A  J for  growth 
along  (1)  to  reach  A  the  change  in  is  also  AJp^  (Tj  and  Aa  are  the  same).  The 
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JpaKbofiow) 

(a) 


AJFB/{bOfiow) 

(b) 

FIG.  12 — Development  of  approximate  normalization  scheme  to  minimize 
dependence  of  growth  constraint  model  on  absolute  values  ofJjc. 
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m  =  dJo/cUpB  ^  =  dJo/dJpB 


JpQ  I  {bOfiow)  '^FB  /  {bOflow) 

(a)  a/W  =  0.1  (b)  a/W  =  0.5 


FIG.  13 — Gradients  (m)  for  stationary  constraint  models  of  SE(B)  specimens 
(a/W  =  0.1,  0.5;  n  =  6,  10;  E/oq  =  500). 


corresponding  changes  in  Jq  values  during  growth  are  shown  as  AJ^  and  Ac/q  .  Us¬ 
ing  the  approximation  that  ai«a2,  Fig.  12  (b)  illustrates  that  the  difference  between 
A  and  A  may  be  given  by  the  expression: 

AJ^  -  AJ®  «  AJpB  •  (mi  -  m2)  . 

Using  this  relationship,  the  constraint  corrected  value  for  point  A  may  be  expressed 
in  terms  of  the  constraint  corrected  value  for  point  B  as  (solve  above  expression  for 
AJp 

=  J™  +  AJ^  “  +  AJf  +  AJpg-  (m,  -  m,)  •  (1*) 

Again, all  J  quantities  in  this  expression  are  understood  to  be  normalized  by  the  ap¬ 
plicable  value  of  bofi  OW 

For  interpretation  of  experimental  data,  the  corrected  toughness  value,  Jq,  is 
sought  given  measured  values  for  Jps  Aa  at  fracture,  in  addition  to  estimated 
values  for  J/c  and  the  tearing  modulus.  To  develop  more  useful  expressions,  it  proves 
convenient  to  define  a  common  reference  position  on  the  stationary  crack  curves.  We 
selected  the  position  where  the  gradient  of  the  stationary  curve  equals  unity. 


mref=  1. 

The  growth  portion  of  each  constraint  model  curve.  At/ o ,  is  normalized  to  define 
AJq  at  the  reference  position  using 

AJq  =  AJq  +  AJ pB  •  (1  ~  m^g)  ^ 

where  m/c  denotes  the  gradient  of  the  stationary  curve  at  J/c  for  the  specific  configu¬ 
ration.  Figures  14  and  15  summarize  results  for  all  the  SE(B)  analyses  (a/W=0.1, 

0.5;  n=5, 10)  for  each  of  the  J/c  and  tearing  modulus  values  using  both  definitions  for 
computation  of  the  growth  crack  curves  (instantaneous  and  cumulative  areas).  Fig¬ 
ure  16  provides  similar  results  for  much  larger  selected  values  of  tearing  modulus. 
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X  10^  ^JjibOfl„)  X  10^ 


(a)  (b) 

tJoKbOfloJ  X  10^  A7o/(fco),„J  X  10’ 


(c) 


(d) 


am  =0.1  ^  =  2.7  ^  =  600 


FIG.  14 — Growth  portion  of  constraint  model  normalized  to  reference  point 
on  stationary  crack  solution  for  SE(B)  (alW=0.1;  n=5,  10). 
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The  normalization  using  Eq  15  accomplishes  the  intended  purpose  of  minimizing  the 
dependence  on  Jjc  for  growth  curves  having  a  common  tearing  modulus.  Our  analy¬ 
ses  also  indicate  a  weak  dependence  of  AJq  on  ai/ao  over  the  range  2.5  ^ai/ao:£  3.1. 

Equations  14  and  15  are  combined  to  yield  the  final  expression  for  Jq  needed  to 
correct  experimentally  measured  toughness  values.  Deleting  the  notation  of  configu¬ 
rations  (1),  (2)  and  the  notion  of  a  specific  point  A  from  Eq  14,  a  more  general  ex¬ 
pression  is  given  by 

Jo  =  +  AJo  . 

(Figure  12  (a)  provides  the  graphical  interpretation  of  this  expression).  Equation  15 
is  solved  for  AJq  and  the  result  substituted  into  the  above  relationship  to  define 

^0  ^0  .  AJq  .  (I  N  (17) 

-  - - j-  — - — - 

^^flow  ^^flow  ^^flow  ^^flow 

where  normalization  of  all  J-quantities  by  the  initial  remaining  ligament 
(6  =  W-  ao)  and  the  flow  stress  is  indicated. 

Although  not  explicitly  used  here,  we  note  the  strong  linear  dependence  of  AJq 
with  AJfb  which  may  lead  to  development  of  useful  expressions  of  the  form 

where  tabular  or  curve  fit  values  for  Jb  can  be  derived  from  Figs.  14  and  15. 

The  following  procedure  summarizes  the  application  of  Eq  17  to  correct  exper¬ 
imental  measures  of  fracture  toughness: 

Quantities  needed: 

■  measured  or  inferred  values  for  n,  Ofiowy  '^J 

■  measured  J  at  cleavage  fracture  (Jfb) 

Steps  in  procedure: 

1.  compute  AJfb  =  JpB  -  Jlc 

2.  look  up  AJq  using  AJfb  iu  the  graphs  provided  in  Figs.  14  and  15 

3.  look  up  m  for  Jfb  =  J/c  using  graphs  in  Fig.  13 

4.  evaluate  using  the  constraint  correction  curves  for  the  stationary  crack 
(  is  Jo  evaluated  at  J/c;  Fig.  1  provides  stationary  crack  solutions.) 

5.  evaluate  Jq  by  substituting  the  above  quantities  into  Eq  17 

Note:  The  experimental  estimates  for  J/c  used  in  this  procedure  should  be  obtained 
by  projecting  early  parts  of  the  i2— curve  back  to  a  suitable  blunting  line,  not  to  a 
0.2mm  offset  line,  to  be  consistent  with  the  finite  element  interpretation  of  Act. 

The  constraint  model  including  crack  growth  expressed  by  Eq  17  and  associated  fig¬ 
ures  retains  the  essential  features  of  original  Dodds-Anderson  model.  This  engineer¬ 
ing  model  captures  the  key  features  of  the  phenomenon  in  a  straightforward  ap¬ 
proach.  Application  of  the  new  model  to  correct  experimental  values  of  fracture 
toughness  for  constraint  involves  only  two  new  steps  compared  to  the  model  without 
crack  growth. 

5.2  Application  of  Constraint  Model  Including  Crack  Growth 

Experimental  data  consisting  of  Jr-Ao  values  [12,17]  that  terminate  in  cleavage 
fracture  are  presented  Figures  17  and  18  for  two  materials,  A515  plate  and  sub- 
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a/lV  =0.5  ^  =  2.7  §-  =  500 


FIG.  15 — Growth  portion  of  constraint  model  normalized  to  reference  point  on 
stationary  crack  solution  for  SE(B)  (alW-0,5;  n=5,  10). 
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FIG.  16 — Growth  portion  of  constraint  scaling  model  normalized  to  reference 
point  on  stationary  crack  solution  for  SE(B)  (alW=0.1;  n=5,  10,  high  tearing 
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merged  arc  welds  in  A533B.  The  A515  data  set  [12]  was  produced  by  testing  SE(B) 
specimens  (a/W  =  0.1)  while  the  A533B  data  set  [17]  was  produced  by  testing  C(T) 
specimens  with  afW  =  0.5.  The  C(T)  specimens  are  analyzed  here  using  computation¬ 
al  results  for  the  deep-notch  SE(B)  specimen.  Table  1  summarizes  the  material 
properties  and  testing  conditions.  Estimates  for  J/c  and  for  the  constant  value  of 
tearing  modulus  indicated  in  the  table  are  found  by  fitting  the  experimental  data 
(note  the  maximum  amount  of  ductile  crack  extension  is  6%  of  bo  for  the  A533B  tests 
and  25%  of  oq  for  the  A515  tests).  The  rightmost  column  in  Table  1  shows  the  ratio 
of  the  tearing  modulus  to  the  the  slope  of  the  blunting  line  expressed  as  a  tearing 
modulus.  Specific  finite  element  analyses  were  performed  for  these  configurations  to 
accommodate  the  actual  tearing  modulus  for  construction  of  the  constraint  model. 


Sources  of  Test 
Results 

ASTM 

Steel 

McCabe  [J7] 

A533B-G 

1 

m 

Kirk,  et  al  [12] 

A515-G70 

20 

5 

296 

545 

88 

420 

0.30 

Table  1-  Values  ofJjc  and  Tj  obtained  by  fitting  experimental  toughness  values. 


Figures  17  (a)  and  18  (a)  compare  measured  toughness  values  expressed  as 
Jp-Aa  curves  with  constraint  corrected  (Jq)  values  obtained  for  each  data  point. 
Filled  squares  denote  measured  toughness  values  (uncorrected).  Filled  triangles  de¬ 
note  e/c -values  corrected  only  for  large-scale  yielding  effects  and  circles  denote 
Jc-values  corrected  for  both  large-scale  yielding  and  crack  growth  effects.  Jo-values 
including  growth  corrections  are  provided  for  the  model  derived  from  the  instanta¬ 
neous  stress  fields  ahead  of  the  growing  crack  and  for  the  model  derived  from  cumu¬ 
lative  areas  within  the  envelope  of  principal  stress  contours.  Jq— values  based  on  the 
cumulative  areas  are  always  larger —  they  include  area  contributions  over  which  the 
principal  stress  equaled  or  exceeded  the  specified  value  during  the  history  of  growth. 

The  correction  for  large-scale  yielding  reduces  the  largest  measured  Jc-values 
by  approximately  a  factor  of  three  for  the  deep— notch  data  set,  see  Fig.  17  (a),  and  by 
a  factor  approaching  five  for  the  shallow-notch  specimens,  see  Fig.  18  (a).  Relative 
to  the  constraint  model  for  large-scale  yielding,  the  new  model  which  includes  both 
growth  and  large-scale  yielding  increases  the  effective  fracture  toughness  (Jq)  by 
approximately  25-30%  for  the  deep-notch  (A533B)  specimens  with  almost  no  in¬ 
crease  due  to  growth  predicted  for  the  shallow-notch  (A515)  specimens.  The  A515 
data  set  has  a  large  tearing  modulus  and  a  high  degree  of  strain  hardening  which 
combine  to  suppress  any  increased  constraint  from  crack  growth.  In  this  data  set, 
Wao/Jic  =  170  and  J/c/ Jo  “2.2;  the  J/?b/ Jo  increases  to  =  5  at  the  last  data  point.  For 
comparison,  Fig.  9  (a)  shows  the  response  of  an  SE(B)  specimen  with  a/W  =  0.1  and 
n  =  5  but  with  a  much  smaller  tearing  modulus  (Tj  =  40  vs.  420  for  the  A515  materi¬ 
al)  and  smaller  deformation  at  the  onset  of  tearing  (Woq/  Jic  =  240  vs.  170  for  the 
A515  material).  In  that  case,  crack  growth  maintains  constraint  constant  at  JpB  I 
«2  under  increased  loading  which  implies  that  during  crack  growth  AJpB  =  A  Jo;  for 
the  A515  data  set,  AJpB  >  A  Jo. 
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Jc  (kJ/nn2)x102  Ln[Ln(^/{1  -  P^))] 


Aa  (mm)  Ln  {Jc)  (kJ/m^) 

(a)  (b) 


A533B  (-15°,0°C),  McCabe  [17] 
C(T),  a/W  =  0.52,  W  =  51  mm,  B=25  mm 

Jlo  =  75  kJ/m2  ^  =  400 


FIG.  17 — Experimental  Jr-^  curves  and  three-parameter  Weibull 
distribution  for  A533B  weld  metal  with  corrections  for  effects  of  crack  growth  and 
large-scale  yielding. 

Figures  17  (b)  and  18  (b)  show  uncorrected  and  corrected  cleavage  toughness  val¬ 
ues  in  the  form  of  Weibull  distributions  for  the  A533B  and  A515  data  sets.  The  Wei¬ 
bull  diagrams  show  three  sets  of  data:  (1)  the  measured  fracture  toughness  values 
iJc),  shown  as  filled  squares;  (2)  toughness  values  corrected  only  for  constraint  (loss) 
due  to  LSY,  shown  as  filled  triangles;  and  (3)  toughness  values  corrected  both  for 
LSY  and  crack  growth  effects,  shown  as  circles.  The  diagrams  illustrate  differences 
in  the  growth  correction  derived  from  instantaneous  and  cumulative  areas  for 
A533B;  for  A515  the  instantaneous  and  cumulative  corrections  are  identical  and 
only  the  cumulative  area  result  is  shown. 

To  construct  the  Weibull  diagrams,  fracture  probabilities  for  the  ranked  tough¬ 
ness  values  are  computed  using  F)  =  (i-0.5)/A(.  In  these  diagrams,  straight  lines  rep¬ 
resent  the  best  fits  to  the  distributions.  The  Weibull  slopes  of  the  dotted  straight 
lines  for  the  uncorrected  distributions  (2.4  and  2.6)  are  slightly  larger  than  the 
theoretical  slope  of  m=2.0  for  toughness  values  expressed  in  terms  of  J.  Similarly, 
slopes  of  the  solid  straight  lines  for  the  corrected  distributions  are  larger  as  well 
(m  =  3.1,  3.6). 

The  trends  of  these  first  application  of  the  growth  correction  model  are  encourag¬ 
ing  with  applications  of  the  model  to  additional  data  sets  required  for  verification. 
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A515(20°C),Kirk,  etal.  [12] 
SE(B),  a/W  =  0.1 ,  W  =  B=51  mm 


Tj  =  420 
n  =  5 


Jic  =  88  kJ/m2 


FIG.  18 — Experimental  Jp-Aa  curves  and  three-parameter  Weibull 
distribution  forA515  with  corrections  for  effects  of  crack  growth  and  large-scale 
yielding. 

especially  data  sets  with  moderate  hardening  and  lower  tearing  modulus  values  rel¬ 
ative  to  the  blunting  line  slope. 

6.  Conclusions 

The  numerical  investigation  of  small  amounts  of  ductile  crack  growth  in  SE(B) 
specimens  relevant  to  cleavage  fracture  described  in  this  work  supports  the  follow¬ 
ing  conclusions: 

1.  The  growing  crack  tip  develops  a  macroscopically  sharp  opening  profile  and  re¬ 
mains  sharp  for  the  maximum  amount  of  growth  considered  in  each  combination 
of  geometry  and  material  properties.  The  crack-tip  opening  angle  remains 
constant  following  the  initial  transient  period  of  growth  in  the  first  1-2  finite 
elements  {Aa^djc),  and  increases  in  proportion  to  increases  in  the  tearing 
modulus. 

2.  During  crack  growth,  J-values  computed  over  domains  very  near  the  tip  ap¬ 
proach  zero.  For  domains  at  distances  from  the  current  tip  greater  than 
0.15-0.20&0,  becomes  path  (domain)  independent  and  agrees  with  J-values 
computed  using  the  procedures  described  in  ASTM  E 1152— 87  (the  deformation 
theory  J  derived  from  growth  corrected  rj  factors). 
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3.  Small  amounts  of  crack  growth  can  impact  significantly  the  stress  and  deforma¬ 
tion  fields  ahead  of  the  extending  tip.  Compared  to  a  stationary  crack  in  an 
SE(B)  loaded  to  the  same  J,  crack  growth  may:  elevate  the  opening  mode 
stresses  on  the  remaining  ligament,  elevate  the  stress  triaxiality  ahead  of  the 
tip,  and  increase  the  absolute  size  of  principal  stress  contours  enclosing  the  tip. 
The  interaction  crack  growth  with  the  degree  of  strain  hardening,  Jic,Tj,a/W 
and  absolute  specimen  size  is  quite  complex.  The  general  effects  of  each  parame¬ 
ter  are;  increased  levels  of  strain  hardening  and  larger  tearing  modulus  values 
suppress  the  effects  of  crack  growth  (drive  the  growth  constraint  model  towards 
the  stationary  crack  model);  large  Jjc  values,  small  a/W  ratios  and  small  speci¬ 
men  sizes  magnify  growth  effects  by  severely  lowering  constraint  at  the  onset  of 
growth. 

4.  The  principal  stress  contours  ahead  of  the  growing  crack  show  a  striking  spatial 
similarity;  while  the  absolute  sizes  of  the  contours  scale  nonlinearly  with  J,  the 
angular  variation  in  shape  undergoes  negligible  changes  during  crack  growth. 
These  same  observations  about  principal  stress  contours  were  made  in  earlier 
work  that  considered  stationary  cracks  in  SE(B)  specimens. 

5.  The  JfbUq  ratios  {FB  denotes  finite-sized  specimens)  predicted  by  the  stress- 
volume  constraint  model  including  crack  growth  show  a  weak  dependence  on 
the  principal  stress  contour  {oiIoq)  selected  for  computation.  For  SE(B)  speci¬ 
mens  with  a/W=0.1,  crack  growth  restores  a  very  weak  dependence  on  o-Joo  rela¬ 
tive  to  the  much  greater  dependence  on  oi/oq  observed  for  a  stationary  crack  at 
high  J/rs-values  (relative  to  specimen  size).  For  SE(B)  specimens  with  a/W=0.5, 
the  dependence  on  oi/oq  develops  for  both  stationary  and  growing  cracks  when 
the  global  bending  field  impinges  too  strongly  on  the  crack  tip  fields. 

6.  A  relatively  simple  normalization  procedure  reduces  the  dependence  of  the 
constraint  model  on  specific  values  of  J/c —  an  important  feature  which  in¬ 
creases  the  utility  of  limited  numerical  solutions.  The  gradient  of  the  response 
(curve),  dJpJdJQ,  for  the  stationary  crack  model  at  J  =  Jjc  provides  sufficient 
information  to  capture  accurately  J/c  effects  on  the  response  during  growth.  De¬ 
pendence  on  the  tearing  modulus  cannot  be  normalized;  however,  the  tearing 
modulus  effects  appear  amenable  to  simple  curve  fitting  of  the  finite  element 
results. 

7.  Cleavage  fracture  values  for  two  different  pressure  vessel  steels  have  been  cor¬ 
rected  for  the  effects  of  ductile  tearing  and  large-scale  yielding  with  the  new 
constraint  model.  The  constraint  model  which  includes  only  large-scale  yielding 
effects  appears  to  overcorrect  measured  Jc— values  in  A533B  that  have  some 
tearing  prior  to  cleavage  fracture  (Jo-values  are  too  small  with  only  the  LSY 
correction).  In  contrast,  shallow  notch  specimens  (afW  =  0.1)  of  a  material 
(A515)  with  high  strain  hardening  (n  =  5)  and  a  large  tearing  modulus 

{Tj  =  400)  exhibit  essentially  no  increase  in  Jo-values  predicted  by  the  growth 
model. 
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